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Boltzmann's principle S(E,N,V) — k\nW relates the entropy to the geomet- 
ric area e s ^ B ' N,v ' of the manifold of constant energy in the (finite- iV)-body phase 
space. From the principle, all thermodynamics and especially all phenomena of 
phase transitions and critical phenomena can be deduced. The topology of the 
curvature matrix C(E, N) (Hessian) of S(E, N) determines regions of pure phases, 
regions of phase separation, and (multi-)critical points and lines. Thus, C(E, N) de- 
scribes all kind of phase-transitions with all their flavor. They are linked to convex 
(upwards bending) intruders of S(E, N), here the canonical ensemble defined by the 
Laplace transform to the intensive variables becomes multi-modal, non-local, and 
violates the basic conservation laws (it mixes widely different conserved quantities) . 
The one-to-one mapping of the Legendre transform gets lost. Within Boltzmann's 
principle, Statistical Mechanics becomes a geometric theory addressing the whole 
ensemble or the manifold of all points in phase space which are consistent with 
the few macroscopic conserved control parameters. Moreover, this interpretation 
leads to a straight derivation of irreversibility and the Second Law of Thermody- 
namics out of the time-reversible microscopic mechanical dynamics. It is the whole 
ensemble that spreads irreversibly over the accessible phase space not the single 
iV-body trajectory. This is all possible without invoking the thermodynamic limit, 
extensivity, or concavity of S(E, N, V) and also without invoking any cosmological 
constraints. Without the thermodynamic limit or at phase-transitions the systems 
are usually not self- averaging, do not have a single peaked distribution in phase 
space. It is further shown that non-extensive Hamiltonian systems at equilibrium 
are described by Boltzmann's principle and not by Tsallis non-extensive statistics. 



PACS numbers: 05.20.Gg,05.70Ln 



I. INTRODUCTION 



There are many attempts to derive Statistical Mechanics from first principles. The earliest are by 
Boltzmann [1-4], Gibbs [5,6], and Einstein [7-10]. The two central issues of Statistical Mechanics 
according to the deep and illuminating article by Lebowitz [11] are to explain how irreversibility 
(the Second Law of Thermodynamics) arises from fully reversible microscopic dynamics, and the 
other astonishing phenomenon of Statistical Mechanics: the occurence of phase transitions. In 
this paper I want to present an easy, straightforward derivation of both aspects directly out of the 
microscopic time-reversal invariant Newton-mechanics invoking a minimum of assumptions. We 
will see how both problems are connected. 

There is an important aspect of Statistical Mechanics which to my opinion was not sufficiently 
considered up to now: Statistical Mechanics and also Thermodynamics are macroscopic theories 
describing the average 1 behavior of all iV-body systems with the same macroscopic constraints. 
It is this fact and nothing else that leads in a simple and straightforward manner to the desired 
understanding of irreversibility, the Second Law for iV-body systems, which obey a completely time 
reversible Hamiltonian dynamics, and leads simultaneously to the full spectrum of phase-transition 
phenomena. It is certainly essential to deduce the Second Law from reversible (here Newtonian) 
and not from dissipative dynamics as is often done because just the derivation of irreversibility from 



1 Here I do not speak of the typical behavior. This would only be the same if the system is self-averaging, 
which I do not demand, see below. 



fully reversible dynamics is the main mystery of Statistical Mechanics. Here a first hint: Whereas 
a single trajectory in the (finite-A^)-body phase space returns after a finite Poincare recurrence 
time a manifold of points develops in general irreversibly with time, see below. 



II. MINIMUM-BIAS DEDUCTION OF STATISTICAL MECHANICS 

Thermodynamics presents an economic but reduced description of a TV-body system with a 
typical size of N ~ 10 23 particles in terms of a very few (M ~ 3 — 8) "macroscopic" degrees 
of freedom (<io/'s). Here we will allow also for much smaller systems of some 100 particles like 
nucleons in a nucleus. However, I assume that always 6N M. The believe that phase transitions 
and the Second Law can exist only in the thermodynamic limit will turn out to be false. 

Evidently, determining only M do/'s leaves the overwhelming number 6N — M dof's undeter- 
mined. All N-body systems with the same macroscopic constraints are simultaneously described 
by Thermodynamics. These systems define an ensemble M of points 2 in the A^-body phase 
space. Thermodynamics can only describe the average behavior of this whole group of systems. 
I.e. it is a statistical or probabilistic theory. Considered on this level we call Thermodynamics 
thermo-statistics or since Gibbs Statistical Mechanics. 

The dynamics of the (eventually interacting) A-body system is ruled by its Hamiltonian H^. 
Let us in the following assume that our system is trapped in an inert rectangular box of volume V 
and there is no further conservation law than the total energy. The motion in time of all points of 
the ensemble follows trajectories in iV-body phase space {qi(t),Pi(t)}\^L 1 (I consider only classical 
mechanics) which will never leave the (6N — l)-dimensional shell (or manifold) £ of constant energy 
E in phase space. We call this manifold the micro- canonical ensemble. 

An important information which contains the whole equilibrium Statistical Mechanics including 
all phase transition phenomena is the area W(E, N) =: e s of this manifold £ in the n-body phase 
space. Boltzmann has shown that S(E, N, V) is the entropy of our system. Thus the entropy and 
with it equilibrium thermodynamics has a geometric interpretation. 

Einstein called Boltzmann's definition of entropy as e.g. written on his famous epitaph 



S=k-lnW 



(1) 



Boltzmann's principle [12] from which Boltzmann was able to deduce thermodynamics. Precisely 
W is the number of micro-states 3 of the A^-body system at given energy E in the spatial volume 
V and further-on I put Boltzmann's constant k = 1: 

W(E,N,V) =tr[e S{E-H N )} (2) 
tr[6(E-H N )} = [ ±(^P) S(E-H N ), (3) 



{qev} Nl V (2^) 3 

eo is a suitable energy constant to make W dimensionlcss, the A" positions q arc restricted to the 
volume V , whereas the momenta p are unrestricted. In what follows, I remain on the level of 
classical mechanics. The only reminders of the underlying quantum mechanics are the measure 
of the phase space in units of 2-kT% and the factor 1/AH which respects the indistinguishability 
of the particles (Gibbs paradoxon). With this definition, eq.(l), the entropy S(E, N,V) is an 
everywhere multiple differentiable, one-valued function of its arguments. This is certainly not the 
least important difference to the conventional canonical definition. 



2 In this paper I denote ensembles or manifold in phase space by calligraphic letters like M. 

3 In the following I will call single points in the 6iV-dim phase-space states or micro-states which are specific 
microscopic realizations of the iV-body system and correspond to single N-body quantum states in quantum 
mechanics. These must be distinguished from macro-states used in phenomenological thermodynamics c.f. 
section V. 
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In contrast to Boltzmann [2,3] who used the principle only for dilute gases and to Schrodinger [13], 
who thought equation (1) is useless otherwise, I take the principle as the fundamental, generic 
definition of entropy. In a recent book [14] cf. also [15,16] I demonstrated that this definition of 
thermo-statistics works well especially also at higher densities and at phase transitions without 
invoking the thermodynamic limit. This is important: Elliot Lieb [17,18] considers the additivity 
of S(E) and Lebowitz [19,11] the thermodynamic limit as essential for the deduction of thermo- 
statistics. However, neither is demanded if one starts from Boltzmann's principle. Boltzmann's 
principle eq.(l) is the only axiomatic assumption necessary for thermo-statistics. 

This is all that Statistical Mechanics demands, no further assumption must be invoked. Neither 
does one need extensivity 4 , nor additivity nor concavity of S(E) cf. [20]. In the next section I 
will show how Boltzmann's principle allows to define phase-transitions in "Small" , non-extensive 
as well in normal "large" extensive systems where our more general definition of phase transitions 
(see below) will coincide with the conventional definition by the Yang-Lee singularities [21,22]. 

Of course one should not wonder if some familiar features of conventional canonical thermo- 
statistics do not hold anymore in "Small" systems. This is discussed in some more detail in 
subsection III A 



III. WHY IS THE MICRO-CANONICAL ENSEMBLE FUNDAMENTAL? 



During the dynamical evolution of a many-body system interacting by short range forces the 
internal energy is conserved. Only perturbations by an external "container" can change the energy. 
I.e. the fluctuations of the energy are 

^ cx y-i/3. (4 ) 

If, however, the diameter of the system is of the order of the range of the force, i.e. the system is 
"Small" , non-extensive, details of the coupling to the container cannot be ignored. 



A. Equivalence of ensembles and self-averaging 

In contrast, the canonical ensemble does not care about these details, assumes the system 
is homogeneous, averages over a Boltzmann-Gibbs (exponential) distribution Pbg{Qc^Poi} = 
_i_^ e -PH{q a ,p a } Q £ ener gy an d fixes only the mean value of the energy by the temperature 1//3. 

In order to agree with the micro, e~^ E W(E) must be sharp in E i.e. self-averaging, what is usu- 
ally not the case in non-extensive systems or at phase transitions of first order. Then one must 
work in the micro, the only orthode ensemble. The micro-ensemble assumes precise - perhaps 
idealized - boundary conditions for each particle independently of whether the system is small 
or large. Therefore, already Gibbs considered the micro-ensemble as the fundamental and the 
canonical as approximation to it. He demonstrates clearly the failure of the canonical in cases of 
phase separation or other situations where both ensemble differ, footnote on page 75 of [5], see 
also [23,24]. 

There are important features where the microcanonical statistics of "Small" systems deviates 
from the "canonical" structure of conventional thermo-statistics of extensive systems in the ther- 
modynamical limit: 

E.g. the familiar Leg endre-trans form structure, a paradigm of "canonical" thermo-statistics, is 
lost. Clearly without self-averaging, fixing an intensive parameter like the temperature T does 



4 Dividing extensive systems into larger pieces, the total energy and entropy are equal to the sum of those 
of the pieces. I will call non-extensive systems where this is not the case in the following also "Small" 
systems [14] (with a capital S\) to stress the paradoxical point that the largest systems in nature (globular 
galaxies) belong to this group as well, nevertheless, they cannot be treated in the thermodynamic limit. 
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not fix the energy sharply. Most evident example is a transition of first order in the canonical 
ensemble at the transition temperature where the energy per particle fluctuates by the specific 
latent heat even in the thermodynamic limit. Related is the occurrence of negative specific heat 
cf. section IV B found in recent experiments on nuclei [25,26] which was predicted many years 
before [27]. Here, there are at least three energies for the same temperature cf. section IV C. 
The present discussions of non-extensive statistics as proposed by Tsallis [28] or recently by Vives 
et al. [29] clearly miss this crucial point. In the Tsallis statistics the entropy is expressed by the 
mean-values of the extensive quantities like <E> [29,30] controlled by a Lagrange parameter j3 or 
(3*. Of course, this is equivalent to the micro-ensemble (is an orthode) only if the variance of the 
energy is small. In one or the other way the thermodynamic limit, extensivity, and self-averaging is 
still demanded where Legendre transformations (may) become one to one. However, in the case of 
non-extensive systems the existence of the thermodynamic limit is unlikely and so is the uniqueness 
of the Legendre transformation. 



B. Tsallis statistics does not apply to Hamiltonian systems at equilibrium 

Tsallis suggested to extend the Boltzmann-Gibbs canonical energy distribution PBG{Qi,Pi} by 
using the q-exponential [31,30]: 

e x q = [l + (l-q)x}^ (5) 

e -f3H{q t , Pi } 

and replacing Pbg = 7^3) ( 6 ) 

N 

where: Z{fi) = I ^ ( ) e"W«> (7) 



J-l 


( d 3 qi d 3 p t 


Nl 


V (2tt^) 3 


-I, 


( d 3 qi d 3 Pl 


Nl 


{ (2wh) 3 



N 

with mean-values <0>bg = I -rW ( ^f' ) 0(q i7 pi) P B G(qi,Pi) (8) 

by: PBctihPi} -» [fUi,Pi}] q 

e ~f3H{q t ,p,} 

where: f{q l ,p l } = q - ^ — (9) 

- lz *^, £ ,i(wj"'*"' (10) 

N 

and calculating mean values <0> q — — - (11) 

D q 

N 

v, = i {qev} — ( ",.>i^r ) o{ qi . Pi )[j( qi . P .]". (12) 



J-l 


( d 3 qi d 3 Pl 


Nl 


\ {2nH) 3 


J-| 


( d 3 qi d 3 p t 


Nl 


V (27rh) 3 


J-| 


( d 3 qi d 3 p l 


Nl 


V (2tt^) 3 



which leads to the Tsallis q-entropy [28,30,32,31]: 



1 _ ^ w Pi 

s = k ^a=i a lim s _ Pa m Pa / 15 n 

q — 1 5->l ^— ' 

The Tsallis g-entropy is very similar to ^-dimension well known in mathematics [33]. Its main 
purpose is to emphasize/suppress small probabilities PJ depending on the parameter q < 1 or 
q > 1 resp.. 

For a closed Hamiltonian system at energy E, the P a are the probabilities for each of the 
W{E) microscopic configurations (quantum states). According to Toral [34] this has of course 



(13) 



p -f3H{q z , Pi } 

Jm[/{«,*}]« - (14) 
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the following consequences: After maximizing S q (E) under variation of P a with the constraint of 
Y^, P a = 1 one obtains the equal probability distribution characterized by Boltzmann's entropy 
W(E) =e s ^: 

e -S(E) ^ = E i_ e (i-9)S(B) 
^otherwise ' S q k q—l ' 

Moreover, following Abe [31] and Toral [34] the original definitions of the microcanonical temper- 
ature and pressure (18) through Boltzmann's entropy S(E, N, V), eq.(l), are the only way within 
Tsallis statistics to define the equilibrium of two systems in weak contact and to fulfill the Zeroth 
Law under energy- and volume exchange see also [29]: 

S(E,N,V) = \nW(E,N,V) (17) 

fdiy 1 dS/dV 
phys ~\dEj phvs ~ dS/dE' 1 ' 

I.e. the physical quantity relevant for equilibrium of Hamiltonian systems, extensive or not, is 
the original Boltzmann entropy S(E) — \n\W(E)\, eq.(l), whatever the non-extensivity index q. 
Therefore, for closed Hamiltonian many-body systems at statistical equilibrium, extensive or not,the 
thermo- statistical behavior is entirely controlled by Boltzmann's principle and the microcanonical 
ensemble as discussed in this paper. Tsallis statistics seems to apply to non-equilibrium situations 
like turbulence [35] or the border of chaos [36] etc. 



IV. PHASE TRANSITIONS WITHIN BOLTZMANN'S PRINCIPLE 

At phase-separation the system becomes inhomogeneous and splits into different regions with 
different structure. This is the main generic effect of phase transitions of first order. Evidently, 
phase transitions are foreign to the (grand-) canonical theory with homogeneous density distri- 
butions. In the conventional Yang-Lee theory phase transitions [21] arc indicated by the zeros 
of the grand-canonical partition sum where the grand-canonical formalism breaks down because 
of the Yang-Lee singularities of the grand-canonical potentials. In the following I show that the 
micro-canonical ensemble gives a much more detailed and more natural insight which moreover 
corresponds to the experimental identification of phase transitions by interfaces (inhomogeneities). 



A. Relation of the topology of S(E, N, V) to the Yang-Lee zeros of Z(T, u, V) 

The grand-canonical partition sum may be obtained out of the micro-canonical one by a double 
Laplace transform. To explore the link to the Yang-Lee singularities I discuss it for the moment in 
the thermodynamic limit (large volume V, large number of particles N but homogeneous constant 
density n = N/V. In this limit it does not matter whether N is discrete or continuous.) 

Z(T, V) = tr — dN e -[E-^N-TS(E)]/T 

JJo eo 

= YL ir de dn e -V[e-»n-Ts { e,n)]/T (jg) 

eo JJo 

v 2 rr°° 

W — // de dn e -^[c°n S t.+lin.+quadr.] 

eo JJo 

and we investigate the specific free energy /(T,/i) = — [ln(Z)]/F in the thermodynamic limit 

V — ► <X>\N/V=const- 

The double Laplace integral (19) can be evaluated asymptotically for large V by expanding the 
exponent as indicated in the third line to second order in Ae, An around the "stationary point" 
e s , n s where the linear terms vanish: 
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1 _ ds 
T ~ ~de 



T 



ds_ 

dn 



(20) 



the only terms remaining to be integrated are the quadratic ones. If the eig en- curvature \\ < 
eq.(24), and eqs.(20) have a single solution (e s ,n s ), this is then a Gaussian integral and yields: 

Z{T,H,V) = — e -V[es-»n s -Ts(e s , ns )]/T // dvi dv2 e V[^vl+\ivl]/l (21) 
e JJ-oo 

= e -F(W)/T (22) 

m „ v) , Sl^Zl _ „ _ _ T „ + IMSMi + o( !^,. (23) 



Here det(e s , n s ) is the determinant of the curvatures (Hessian) of s(e, n), v±, t>2 are the eigenvectors 
of the Hessian. 



= AiA 2 , Ai > A 2 (24) 



Ai can be positive or negative. If Ai < and eqs.(20) have no other solution, the last two terms 
in eq.(23) go to 0, and we obtain in the thermodynamic limit (V — > oo) the familiar result for the 
free energy density: 

f(T, (j,, V — ► oo) = e s - im s - Ts s . (25) 

I.e. the curvature Ai of the entropy surface s(e, n, V) or the largest eigenvalue of the Hessian matrix 
decides whether the grand- canonical ensemble agrees with the fundamental micro- ensemble in the 
thermodynamic limit. If this is the case and eqs.(20) have a single solution or s(e,n) touches its 
concave hull at e s ,n s , then there is a pointwise one to one mapping of the micro-canonical entropy 
s(e,n) to the grand-canonical partition sum Z(T,fI), and hi[Z(T, jJL)]/V or f(T, /x) is analytical 
in z = e^. Due to Yang and Lee we have then a single, stable phase [22]. Otherwise, the 
Yang-Lee zeros of Z(T,fi) reflect anomalous points/regions of Ai > (det(e,n) < 0) where the 
canonical partition sum does not reflect local properties of the micro-ensemble, i.e. does not respect 
the conservation laws, and mixes conserved quantities. This is crucial. As det(e s , n s ) can be studied 
for finite or even small systems as well, this is the only proper extension of phase transitions to 
"Small" systems. 



det(e, n) 



d 2 s 



d s d s 
dedn dn 2 



B. The physical origin of the wrong curvature 

I will now discuss the physical origin of the convex (upwards bending) intruders in the entropy 
surface for systems with short-range coupling in two examples. 
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1. Liquid-gas transition in sodium clusters 



course with fluctuations. The latent heat per atom is 
qiat = ez — e\. Attention: the curve s(e) is artifically 
sheared by subtracting a linear function 25 + e * 11.5 in 
order to make the convex intruder visible. s(e) is al- 
ways a steeply monotonic rising function.We clearly see 
the global concave (downwards bending) nature of s(e) 
and its convex intruder. Its depth is the entropy loss due 
to the additional correlations by the interfaces. It scales 
oc iV -1 / 3 . From this one can calculate the surface ten- 
sion per surface atom a sur f/Tt r = As sur f * No/N sur f. 
The double tangent (Gibbs construction) is the concave 
hull of s(e). Its derivative gives the Maxwell line in the 
caloric curve T(e) at T tr - In the thermodynamic limit the 
intruder would disappear and s(e) would approach the 
double tangent from below. Even though, however, the 
probability of configurations with phase-separations are 
suppressed by the (infinitely small) factor e~ N /3 relative 
to the pure phases and the distribution remains strictly 
bimodal in the canonical ensemble. 

In tabic (I) I compare the "liquid-gas" phase transition in sodium clusters of a few hundred 
atoms with that of the bulk at 1 atm. c.f. also fig.(l). 

TABLE I. Parameters of the liquid-gas transition 
of small sodium clusters (MMMC-calculation) in com- 
parison with the bulk for rising number No of atoms, 
N aU rf is the average number of surface atoms (es- 
timated here as J2 ^cluster) °f au clusters together. 
o /T tr = As sur f * No/N sur f corresponds to the surface 
tension. Its bulk value is adjusted to the experimental 
input values used for the binding energies of clusters as 
given by Brechignac et al. c.f. [14]. 

Conclusion: For systems with short range interactions a convex intruder in s(e, n) appears with 
the fragmentation of the system into several clusters and monomers. The number of surface par- 
ticles scales with the depth of the intruder (surface entropy). I.e. the convex intruder signals the 
preference of the system to become inhomogeneous, the characteristic signal for the separation of 
different phases (liquid and gas) at a phase transition of first order. 

2. The topology of the entropy surface S(E, N) for Potts lattice gases 

Having discussed in the previous example a system with a single thermodynamic degree of 
freedom or control parameter (the energy E) we will now study more subtle features. If the system 
has two, or more, degrees of freedom, e.g. energy E = Ve and particle number N = Vn, where 
V is the volume, we can have phase boundaries and critical points. This is similar to the classical 
P — V diagram of the liquid-gas phase transition in the Van der Waals theory. We are now able to 
identify multi-critical points. These were previously studied in the canonical ensemble only, where 
sophisticated finite size scaling is needed to identify these points. 

As example we investigate the 3-states diluted Potts model on a finite 2-dim (here L 2 — 50 2 ) 
lattice with periodic boundaries in order to minimize effects of the external surfaces of the system. 
The model is defined by the Hamiltonian: 




0.3 0.5 0.7 0.9 1.1 1.3 



FIG. 1. MMMC [14] simulation of the 
entropy s(e) per atom (e in eV per atom) 
of a system of No = 1000 sodium atoms 
with an external pressure of 1 atm. At the 
energy ei the system is in the pure liquid 
phase and at e^ in the pure gas phase, of 



iVo 



200 



1000 



3000 



bulk 



Na 



Ttr [K] 



940 



990 



1095 



qiat [eV] 



0.82 



0.91 



0.94 



Sbo t l 



10.1 



10.7 



9.9 



As sur 



0.! 



0.56 



0.44 



N 



surf 



39.94 



98.53 



186.6 



Q~JTu 



2.75 



5.68 



7.07 



1156 



0.923 



9.267 



7.41 
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n.n. pairs 

H = - (26) 
n = L~ 2 N = L- 2 ^2ot. 

i 

Each lattice site i is either occupied by a particle with spin Cj = —1,0, or 1, or it is empty 
(vacancy). The sum is over pairs of neighboring lattice sites and the occupation numbers are: 

J 1 , spin particle in site i , 07 -. 
1 \ , vacancy in site i 

This model is an extension of the ordinary (q — 3)-Potts model to allow also for vacancies. At 
zero concentration of vacancies (n — 1), the system has in the thermodynamic limit a continuous 
phase transition at e c = 1 + -i= w 1.58 [37,38]. With rising number of vacancies the probability to 
find a pair of particles at neighboring sites with the same spin orientation decreases. I.e. this is 
similar to a larger number of spin orientation q e f t on each lattice site in the ordinary Potts model, 
where we know that the transition of second order becomes a transition of first order for q > 4. 
Thus the inclusion of vacancies has the effect of an increasing effective q e f / > 3. This results in an 
increase of the critical energy of the continuous phase transition with decreasing n and provides 
a line of continuous transition, which is supposed to terminate when q e ff becomes larger than 4. 
From here on the transition becomes first order. 

At smaller energies the system is in one of three ordered phases (spins predominantly parallel 
in one of the three possible directions). 

Figure (2) shows how for a small system of 50*50 lattice points all phenomena of phase transitions 
can be studied from the topology of the determinant of curvatures (24) in the micro-canonical 
ensemble. In this example the second curvature A2 < so sign(det) = — sign(Ai). 




-2 -1.5 -1 -0.5 

e 

FIG. 2. Conture plot of the curvature determinant, eq.(24), of the 2-dim Potts-3 lattice gas with 50 * 50 
lattice points, n is the number of particles per lattice point, e is the total energy per lattice point; Dark 
grey lines: det = 0, boundary of the region of phase coexistence (det < 0), in the triangle AP m B; Light 
grey line: minimum of det(e, n) in the direction of the largest curvature (vA max ■ V det = 0), line of second 
order transition; In the triangle AP m C pure ordered (solid) phase (det > 0); Above and right of the line 
CP m B pure disordered (gas) phase (det > 0); The crossing P m of the boundary lines is a multi-critical 
point. It is also the critical end-point of the region of phase separation (det < 0). The light gray region 
around the multi-critical point P m corresponds to a flat region of det(e,n) ~ i.e. V Ai~ 0, details see 
[15]; C is the critical point of the ordinary Potts, q — 3 model (without vacancies) in the thermodynamic 
limit. 



C. Systematics of phase transitions in the micro-ensemble 

Now we can give a systematic and generic classification of phase transitions in terms of the 
topology of curvatures of s(e,n) which applies also to "Small" systems: 



S 




order j 


arameter 















order parameter 



A single stable phase by Ai < 0. Here s(e,n) is lo- 
cally concave (downwards bended) in both directions 
and eqs.(20) have a single solution e s ,n s . Then there 
is a one to one mapping of the grand-canonical <->the 
micro-ensemble. The order parameter is the direction 
vi of the eigenvector of largest curvature Ai. In many 
situations one may have only locally Ai < 0, however, 
there may be further solutions to eqs.(20) farther away. 
In such cases which have no equivalent in the canonical 
ensemble, we will still speak of regions in {e,n} of pure 
phases embedded in regions of phase- separation. 



A transition of first order with phase separation and 
surface tension (c.f. subsection IV B 1) indicated by Ai > 
0. s(e, n) has a convex intruder (upwards bended) in the 
direction vi of the largest curvature. Then eqs.(20) have 
multiple solutions, at least three. The whole convex 
area of {e,n} is mapped into a single point (T, /i) in 
the grand-canonical ensemble (non-locality). I.e. if the 
largest curvature of S(E, N) is Ai > both ensembles 
are not equivalent, the (grand-) canonical ensemble is 
non-local in the order parameter and violates basic 

conservation laws. See also [14-16,39]. The region in the plane of conserved control- 
parameters e, n where we have a separation of different phases, Ai(e,n) > 0, is bounded 
by lines with Ai(e, n) — 0. On this boundary is the end-point of the transition of first order, 

• where we have a continuous ("second order") transition with vanishing surface tension, 
where two neighboring phases become indistinguishable. This is at points where the two 
stationary points move into one another to become the critical end-point of the first order 
transition. This is then also a maximum of Ai. I.e. where Ai(e, n) = and v> 1=0 • VAi = 0. 
These are the catastrophes of the Laplace transform E — > T. Here va 1= o is the eigenvector 
of the Hessian belonging to the largest curvature eigenvalue Ai = 0. (vi plays the role of 
the order parameter of the transition. In this direction one moves fastest from one phase 
to the other.) At the other points of Ai(e,n) = one of the two coexisting phases gets 
depleted. Furthermore, there may be also whole lines of second-order transitions like in the 
anti-fcrro-magnetic Ising model c.f. [14]. 

• Finally, a multi-critical point where more than two phases become indistinguishable is at 
the branching of several lines in the {e, n}-phase-diagram with Ai = 0,V Ai= 0. Fig. 2 gives 
an illustration of a multicritical point in a small system. 

Just a comment on the interesting suggestion by Pettini et al. [40] to characterize phase tran- 
sitions by topological changes of the potential energy: They claim that phase transitions are 
signalized by sudden changes of the topology of the configurational energy. 

However, one has to keep in mind: In the thermodynamic limit, first order transitions are 
characterized by two essential things: 

1. There are two or more phases (different types of configurations, liquid and gas, or ordered 
and disordered) which coexist at the same temperature T = [dS/dE] . 

2. Between both configurations the total energy per particle, or the potential energy per particle 
differs by the specific latent heat. 

Of course not every topological jump in the configuration energy corresponds to a phase transi- 
tion [40] and not every jump is accompanied by the same kinetic energy per particle (temperature) 
or more precisely by the same slope T — [dS/dE]^ 1 of the entropy. Evidently, topological changes 
are likely neccessary conditions (Pettini) for phase transitions to occur, however, the sufficient 
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condition is the same temperature T — [dS/dE] -1 at two different values of the order-parameters, 
energy or magnetization etc. 

Also the suggestion of Borrmann et al. [41] to characterize phase transitions in small systems by 
the zeros of the canonical partition sum must be commented: This idea extends the idea of Yang 
and Lee and the work of Grossmann [42]. First, the canonical partition sum is a derived quantity 
which is not an orthode at phase transitions, c.f. section III, as was already noticed by Gibbs, 
footnote on page 75 of [5]. Moreover, in some situations of non-extensive systems it may not even 
exist. 



V. MEASURING A MACROSCOPIC OBSERVABLE 

After succeeding to deduce equilibrium statistics including all phenomena of phase transitions 
from Boltzmann's principle alone even for "Small" systems, i.e. non-extensive many-body sys- 
tems [14], it is challenging to explore how far this "most conservative and restrictive way to ther- 
modynamics" [43] is able to describe also the approach of (possible "Small" ) systems to equilibrium 
and the Second Law of Thermodynamics. 

Thermodynamics describes the development of macroscopic features of many-body systems with- 
out specifying them microscopically in all details. Before I address the Second Law, I have to clarify 
what I mean with the label "macroscopic observable" . 

A single point {Qi(t),Pi(t)}i—i t ... t N in the A-body phase space corresponds to a detailed speci- 
fication of the system with all degrees of freedom (do/'s) completely fixed at time t (microscopic 
determination, the curly brackets indicate the whole set of 6 A coordinates q%,p% of all particles i). 
Fixing only the total energy E of an A-body system leaves the other (6 A — l)-dcgrees of freedom 
unspecified. A second system with the same energy is most likely not in the same microscopic 
state as the first, it will be at another point in phase space, the other dof's will be different. I.e. 
the measurement of the total energy Hn, or any other macroscopic observable M, determines a 
(6A — l)-dimensional sub-manifold £ or M in phase space. (The manifold M is called by Lebowitz 
a macro-state [19,11] which contains Tm — W(M) micro-states. I, however, prefer to use the name 
"state" only for micro-states or points in phase space.) All points (the micro-states) in A-body 
phase space consistent with the given value of E and volume V, i.e. all points in the (6A — 1)- 
dimcnsional sub-manifold £(N, V) of phase space are equally consistent with this measurement. 
£(N, V) is the micro-canonical ensemble. This example tells us that any macroscopic measurement 
is incomplete and defines a sub-manifold of points in phase space not a single point. An additional 
measurement of another macroscopic quantity B{q,p} reduces £ further to the cross-section £ n B, 
a (6A — 2)-dimensional subset of points in £ with the volume: 

W(B, E,N,V) = ±J {^0) toS(E - H N {q,p}) S(B - B{q lP }) (28) 

If H^{q,p\ as well as also B{q,p\ are continuous diffcrcntiablc functions of their arguments, what 
I assume in the following, £ n B is closed. In the following I use W for the Riemann or Liouvillc 
volume (Hausdorff measure) of a many- fold. 

Microcanonical thermostatics gives the probability P(B, E, A, V) to find the A-body system in 
the sub-manifold £ n B(E, A, V): 

P(B E N V)- W ( B > E > N > V ) _ HW(B,E,N,V)]-S(E,N,V) (2q] 

( ' ' ' ] ~ W(E,N,V) [ ' 

This is what Krylov seems to have had in mind [44] and what I will call the "ensemble probabilistic 
formulation of Statistical Mechanics (EPS) " . 

Similarly thermodynamics describes the development of some macroscopic observable B{q t ,pt} 
in time of systems which were specified at an earlier time to by another macroscopic measurement 
A{qo,po}- It is related to the volume of the sub-manifold M(t, to) = A(to) n B(t) n £: 

W(A,B,E,t) = ^ J ( f^Kp ) S ( B -B{quPt})S(A-A{q ,po})e a 6(E-H{q u p t }), 

' (30) 
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where {qt{qo,Po},Pt{qo,Po}} is the set of trajectories solving the Hamilton-Jacobi equations 



dH 



dH 



i = 1---N 



(31) 




with the initial conditions {q(t = t ) = q n ; p(t = t ) = Po}- For a large system with N ~ 10 23 
the probability to find a given value B(t), P(B(tj), is usually sharply peaked as function of B at 
its typical value. Such systems are called self-averaging. Ordinary thermodynamics treats systems 
in the thermodynamic limit N — > oo and gives only <B(t)>. However, here we are interested 
to formulate the Second Law for "Small" systems i.e. we are interested in the whole distribution 
P(B(t)) not only in its mean value <B(t)>. There are also many situation where the system 
is not self-averaging, where a finite variance remains even in the thermodynamic limit. (E.g. at 
phase transitions of first order the energy per particle fluctuates in the canonical ensemble by the 
specific latent heat.) 

There is an important property of macroscopic measurements: Whereas at finite times Hamilton 
dynamics evolves a compact region of phase space again into a compact region, this does not need 
to be so at infinite times. Then, at t — > oo, the set may not be closed anymore (perhaps a fractal, 
see below). This means there exist series of points {a n } G A(t = oo) which converge to a point 
lim„_ 

oo ei n — : a n— oo which is not in Ait — oo). E.g. such points ci n — oo ^ *4(oo) may have intruded 
from the phase space complementary to A(t ). Illustrative examples for this evolution of an initially 
compact sub-manifold into a fractal set are the generalized baker transformations discussed in this 
context by ref. [45,46]. See reference [33] for the fractal distribution produced by the general 
baker transformation. (As any housewife knows, a baker dough becomes an infinitely thin (fractal) 
puff pastry after pounding and folding it infinitely often.) Then no macroscopic (incomplete) 
measurement can resolve a n=00 ^ A(t = oo) from its immediate neighbors a n G A(t = oo) in 
phase space with distances \a n — a n=co \ less then any arbitrary small S. In other words, at the 
time t — to — > oo no macroscopic measurement with its incomplete information about {<Zt=oo,f>t=oo} 
can decide whether {qo{qt=oo,Pt=oo},Po{qt=oo,Pt=oo}} G A(to) or not. I.e. any macroscopic theory 
like thermodynamics can only deal with the closure of A(t — ► oo). If necessary, the sub-manifold 
A(t — > oo) must be artificially closed 5 to A(t = oo) as developed further in section VII. Clearly, 
in this approach this is the physical origin of irreversibility. 

VI. ON EINSTEIN'S OBJECTIONS AGAINST THE EPS-PROBABILITY 

Before I proceed I must comment on Einstein's attitude to the principle [47]: Originally, Boltz- 
mann called W the "Wahrscheinlichkeit" (probability), i.e. the relative time a system spends (along 
a time-dependent path) in a given region of QN-dim. phase space. Our interpretation of W to be the 
number of "complexions" (Boltzmann's second interpretation) or quantum states (trace) with the 
same energy was criticized by Einstein [12] as artificial. It is exactly that criticized interpretation 
of W which I use here and which works so excellently [14]. 

According to Abraham Pais: "Subtle is the Lord" [47], Einstein was critical with regard to 
the definition of relative probabilities by eq.(29), Boltzmann's counting of "complexions". He 
considered it as artificial and not corresponding to the immediate picture of probability used in the 
actual problem: "The word probability is used in a sense that does not conform to its definition as 
given in the theory of probability. In particular, cases of equal probability are often hypothetically 
defined in instances where the theoretical pictures used are sufficiently definite to give a deduction 
rather than a hypothetical assertion" [12]. He preferred to define probability by the relative time 
a system (a trajectory of a single point moving with time in the V-body phase space) spends in 
a given subset of the phase space. However, is this really the immediate picture of probability 
used in Statistical Mechanics? This definition demands the ergodicity of the trajectory in phase 
space. As I discussed above, thermodynamics as any other macroscopic theory handles incomplete, 



5 First t — > oo then the closure, not the other way round c.f. however, the discussion in VIII. 
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macroscopic informations of the iV-body system. It handles, consequently, the temporal evolution 
of finite sized sub-manifolds - ensembles - not single points in phase space. In the case of a very 
large system the typical outcomes of macroscopic measurements are calculated. Nobody waits in 
a macroscopic measurement, e.g. of the temperature of a gas, long enough that an atom can cross 
the whole system. 

In this respect, I think the EPS version of Statistical Mechanics is closer to the experimental 
situation than the duration-time of a single trajectory. Moreover, in an experiment on a small 
system like an excited nucleus, which then may fragment statistically later on, the average over a 
multiple repetition of scattering events is taken. No ergodic covering of the whole phase space by a 
single trajectory in time is demanded. Fragmenting nuclei at such high excitation have a too short 
lifetime. This is analogous to the statistics of a falling ball on a Galton's nail-board where also a 
single trajectory does not touch all nails but is random. Only after many repetitions the smooth 
binomial distribution is established. For the discussion of the Second Law in finite systems, this is 
the correct scenario, not the time average over a single ergodic trajectory. 



VII. FRACTAL DISTRIBUTIONS IN PHASE SPACE, SECOND LAW 

Here I will first describe a simple working-scheme (i.e. a sufficient method) which allows to 
deduce mathematically the Second Law. Later, I will show how this method is necessarily implied 
by the reduced information obtainable by macroscopic measurements or theories. 

Let us examine the following Gedanken experiment: Suppose the probability to find our system 
at points {qtTPt]\ in phase space is uniformly distributed for times t < to over the sub-manifold 
£(N,Vi) of the TV-body phase space at energy E and spatial volume V\. At time t > to we 
allow the system to spread over the larger volume V2 > V\ without changing its energy. If the 
system is dynamically mixing, the majority of trajectories {qt,Pt}^ in phase space starting from 
points {qo,Po}i with qo <E V\ at t will now spread over the larger volume V2. Of course the 
Liouvillean measure of the distribution A4(t, to) in phase space at t > to will remain the same 
(= tr{£(N, Vi)]) [48]. (The label {qo G V\} of the integral means that the positions {qo}^ are 
restricted to the volume V\, whereas the momenta {po}\ are unrestricted.) 

1. (33) 

But as already argued by Gibbs the distribution M.(t, to) will be filamcntcd like ink in water and 
will approach any point of £ (N, V2) arbitrarily close, limt_>oo A4(t, to) becomes dense in the new, 
larger £(N, V2). The closure A4(t = 00, to) becomes equal to £(N, V2). This is clearly expressed 
by Lebowitz [19,11]. 

In order to express this fact mathematically, I transform integrals over the phase space like (3) 
or (28): 

W(E,N,t,t ) = ±-f (lOTn) e 6(E-H N {q uPt }) (34) 

into: 



tr[M(t,t )}\ {qoeVl} = 



because of: 

o{qo,Po\ 
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J (d 3 q t d 3 p t ) N ■ ■ ■ = J da!--- da 6N ■ 



da 6 N := 



E 



dH 



dH 



dq t dpi 



l|Vff|| 
W(E,N,t,t ) 



E 



<9# 



E 



dH 

dpt , 



N\{2mi r - J{ qo { quPt } CVl } 
Now, I redefine Boltzmann's definition of entropy eq.(l to 3): 



du\- ■ ■ do-QN-i 



l|Vff|| 



S = ]n(W(E, N, V) 
W(E,N,V)=tr[e 6(E-H N )] 



tr[8(E - H 



N 



by replacing the Riemannian integral for W by its box-counting "measure" : 



W(E,N, V) 



1 



d 3 q t d 3 p t 

(27T?i) 3 



N 



e S(E-H N ), 



(35) 
(36) 

(37) 
(38) 



(39) 
(40) 

(41) 



(42) 



{qo{qt,pt}eVi} 

i.e. the volume of M. by that of its closure Ai. In detail we perform the following steps: 

M 5 (t,t ) :=<G> S *vol box j[M{t,t )}, (43) 

to obtain volbox,s[-M.(t, to)] we cover the d-dim. sub-manifold A4(t,to), here with d — (67V — 1), of 
the phase space by a grid with spacing 5 and count the number N$ oc S~ d of boxes of size S 6N , which 
contain points of M{t,t ). This is illustrated by fig. (3). Then vol box ^[M(t,t )] := 6 d Ng[M(t, to)] 



and <G>s is the average of 



over these non-empty boxes of size S. 



N\(2wh)3 N \\VH\\ 

The liniA^nvoUnf ,s[A4(t, tp)] is the box-counting volume of M(t,to) which is the same as the 



volume of its closure M(t, t ), see below: 

vokoz[M(t,to)] ~]ms^oS d Ns[M{t,to)] 

with lim * = inf [lim *] and write symbolically: 

KmM s (t,t )=:£ ± ff^Pi\ N £oS{ E-H 



(44) 
(45) 



where ijd means that this integral should be evaluated via the box-counting volume (the limit of 
expression 43 with the use of 44) here with d — 6N — 1. This is illustrated by the following figure. 
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V a V b ► Va + Vb 




t<to t>t 
FIG. 3. The compact set M(to), left side, develops into an increasingly folded "spaghetti" -like distri- 
bution M(t, to) in the phase-space with rising time t. The right figure shows only the early form of the 
distribution. At much later times it will become more and more fractal and finally dense in the new phase 
space. The grid illustrates the boxes of the box-counting method. All boxes which overlap with M(t,to) 
contribute to the box-counting volume vol(, 0:Ei( 5 and are shaded gray. Their number is Ng 



The volume of phase space covered by M$(t, to) is > W(E,N,Vi). For finite times because of 
Liouville's theorem eq.(33), see also section (V), we have 

Me-+oM s (t,to) = W(E,N,t ,t ) = W{E,N,V X ) (46) 

At t — > oo the two limits 5 — > 0, t — > oo do in general not commute and as assumed by Gibbs the 
manifold M(t — > oo) becomes dense in the new micro-canonical manifold Then 

Hm^o lim M s(t, to) = W{E, N, V 2 ) > W{E, N, V x ). (47) 

t— >oo 

This is the Second Law of Thermodynamics. The box-counting is also used in the definition of the 
Kolmogorov entropy, the average rate of entropy gain [49,33]. 

The box-counting "measure" is analogous to the standard method to determine the fractal 
dimension of a set of points [49] by the box-counting dimension: 

dinW.M(Mo)] :=l^ J nNs[Mit / 0)] . (48) 

— mo 

Like the box-counting dimension, the box-counting "measure" has the peculiarity that it is equal 
to the measure of the smallest closed covering set. E.g.: The box-counting volume of the set of 
rational numbers volb ox {Q} between and 1 is = 1, and thus equal to the measure of the real 
numbers, c.f. Falconer [49] section 3.1. This is the reason why the box-counting "measure" is not 
a measure in its mathematical definition because then wc should have 



vol 



box 



E cm 

iC{Q} 



^ vok ox [M i ] = 0, (49) 

iC{Q} 



therefore the quotation marks for the box-counting "measure" , c.f. appendix X. 

Coming back to the the end of section (V), the volume W(A, B, ■ ■ ■ , t) of the relevant ensemble, 
the closure M(t,to) must be "measured" by something like the box-counting "measure" (42) with 

the box-counting integral ijd, which must replace the integral in eq.(3). Due to the fact that the 

box-counting volume is equal to the volume of the smallest closed covering set, the new, extended, 
definition of the phase-space integral eq.(42) is for compact sets like the equilibrium distribution 
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£ identical to the old one eq.(3) and our redefinition of the phase-space integral by box-counting 
changes nothing for equilibrium statistics. Therefore, one can simply replace the old Boltzmann- 
dcfinition of the number of complexions and with it of the entropy by the new one (42) of course 
with the understanding that the closure operation with the box-counting volume (42) should be 
done after the times were specified. 

VIII. CONCLUSION 

In this paper I showed that Boltzmann's principle eq.(l) covers in a simple and straight way 
both of Lebowitz's central issues of statistical mechanics [11]. Earlier formulations of these ideas 
can be found in [50,51]. Lebowitz emphasises the neccessity of self-averaging for thermodynamics 
which describes the typical outcome of a macroscopic measurement. This can only be expected for 
large systems, in the thermodynamic limit. However, there are many situations where even large 
systems are not sef- averaging. E.g. at phase transitions of first order. Moreover, a whole world of 
non-extensive systems, like the "Small" systems, show broad, often not single peaked, phase-space 
distributions. E.g. in scattering experiments on nuclei or atomic clusters an average over millions 
of events is taken. Thus the whole distribution in the accessible phase space is measured. These 
are certainly the most interesting situations. An extension of statistical mechanics to cover also 
these is demanded. 

Macroscopic measurements M determine only a very few of all 67V dof . Any macroscopic theory 
like thermodynamics deals with the area M of the corresponding closed sub-manifolds M. in the 
6iV-dim. phase space not with single points. The averaging over ensembles or finite sub-manifolds 
in phase space becomes especially important for the micro-canonical ensemble of a finite or any 
other not self-averaging system. 

Because of this necessarily coarsed information, macroscopic measurements, and with it also 
macroscopic theories are unable to distinguish fractal sets M from their closures M. Therefore, I 
make the conjecture: the proper manifolds determined by a macroscopic theory like thermodynam- 
ics are the closed M. However, an initially closed subset of points at time to does not necessarily 
evolve again into a closed subset at t = oo and the closure operation must be explicitely done after 
setting the times in order to obtain a quantity that is relevant for a macroscopic theory and can be 
compared to thermodynamics. As the closure operation and the t — > oo limit do not commute, the 
macroscopic dynamics becomes irreversible. 

Here is the origin of the misunderstanding by the famous reversibility paradoxes which were 
invented by Loschmidt [52] and Zcrmclo [53,54] and which bothered Boltzmann so much [55,56]. 
These paradoxes address to trajectories of single points in the iV-body phase space which must 
return after Poincare's recurrence time or which must run backwards if all momenta are exactly 
reversed. Therefore, Loschmidt and Zermclo concluded that the entropy should decrease as well 
as it was increasing before. The specification of a single point in QN-dim phase-space and the 
reversion of all its 3N momentum components demands of course a microscopic exact specification 
of all 6N degrees of freedom not a determination of a few macroscopic degrees of freedom only. 
No entropy is defined for a single point. Thermodynamics is addressed to the whole manifold, 
ensemble of systems with the same macroscopic constraints. The ensemble develops irreversibly 
even though the underlining Newtonian dynamics of each phase-space point is fully reversible. It 
is highly unlikely that all points in the ensemble A4(t, to) have commensurable recurrence times so 
that they can return simultaneously to their initial positions. Once the manifold has spread over 
the larger phase space it will never return. 

Also other misinterpretation of Statistical Mechanics are pointed out: The existence of phase 
transitions and critical phenomena are not linked to the thermodynamic limit. They exist clearly 
and sharply in "Small" , non-extensive systems as well. Closed non-extensive Hamiltonian systems 
at equilibrium do not follow Tsallis non-extensive statistics [28]. Boltzmann's principle describes 
the equilibrium and the approach of the equilibrium of extensive as well of non-extensive Hamil- 
tonian systems. 

By our derivation of micro-canonical Statistical Mechanics for finite, eventually "Small" systems 
various non-trivial limiting processes are avoided. Neither does one invoke the thermodynamic limit 
of a homogeneous system with infinitely many particles nor does one rely on the ergodic hypothesis 
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of the equivalence of (very long) time averages and ensemble averages. As Bricmont [43] remarked 
Boltzmann's principle is the most conservative way to Thermodynamics but more than that it is 
the most straight one also. The single axiomatic assumption of Boltzmann's principle, which has a 
simple geometric interpretation, leads to the full spectrum of equilibrium thermodynamics including 
all kinds of phase transitions and including the Second Law of Thermodynamics. 

In this paper, I take the fact serious that Thermodynamics as well as any other macroscopic 
theory handles ensembles or sub-manyfolds and not single points in phase-space. Thus the use of 
ensemble averages is justified directly by the very nature of macroscopic (incomplete) measurements. 
Entropy s(e, n) is the natural measure of the geometric size of the ensemble. With the Boltzmann 
definition of s(e, n) Statistical Mechanics becomes a geometric theory. The topology of its curvature 
indicates all phenomena of phase transitions independently of whether the system is "Small" or 
large. 

Coarse-graining appears as natural consequence of the ensemble-nature. The box-counting 
method mirrors the averaging over the overwhelming number of non-determined degrees of free- 
dom. Of course, a fully consistent theory must use this averaging explicitly. Then one would not 
depend on the order of the limits lfrn^o lim^oo as it was assumed here. Presumably, the rise of 
the entropy can then already be seen at finite times when the fractality of the distribution in phase 
space is not yet fully developed. The coarse-graining is no more a mathematical ad hoc assump- 
tion. It is the necessary consequence of the averaging over the 6N — M uncontrolled degrees of 
freedom. Moreover the Second Law in the EPS-formulation of Statistical Mechanics is not linked 
to the thermodynamic limit as was thought up to now [11,19]. 

In this paper I did not contribute anything to the problem of describing irreversible thermo- 
dynamics of stationary dissipative systems as it is discussed e.g. by Gilbert and Dorfman [57,46], 
Rondoni and Cohen [58]. As mentioned already dissipation does not exist in the microscopic dy- 
namics. It is not clear to me how far the inclusion of dissipation predefines the arrow of time 
already which should have been deduced from the theory. The main problem for me was the 
derivation of irreversibility from fully time reversible microscopic dynamics under maximally clear 
conditions, i.e. of a microcanonical closed system. 

Gaspard [59,60] considers systems obeying a dynamics that preserves the phase-space volume, 
i.e satisfying Liouville's theorem, but under non-equilibrium steady state conditions. Similarly to 
the present approach he had to coarse grain (width S) the accessible phase space. In conformity 
to the standard view of thermodynamics being based on the thermodynamic limit [19] he then 
proves the rise of the entropy after the limits (in that order) 5 — ► 0, V — > oo. However, in this 
limit also the Poincare recurrence time becomes infinite and Zcrmclo's piercing argument becomes 
blunted. So in this approach Gaspard cannot treat our problem of the Second Law in a finite 
closed Hamiltonian system which seems to me to be the heart of the reversibility paradox. 

IX. ACKNOWLEDGMENT 

Thanks to A.Ecker for mathematical advises and to J.Barrc, J.E.Votyakov and V.Latora who 
had useful suggestions to improve the text. I am grateful to S.Abe, M.Baranger, E.G. D.Cohen, 
P.Gaspard, R.Klages, A.K. Rajagopal and M.Pettini for many illuminating discussions. 

X. APPENDIX 

In the mathematical theory of fractals [49] one usually uses the Hausdorff measure or the Haus- 
dorff dimension of the fractal [33]. This, however, would be wrong in Statistical Mechanics. Here 
I want to point out the difference between the box-counting "measure" and the proper Hausdorff 
measure of a manifold of points in phase space. Without going into too much mathematical details 
I can make this clear again with the same example as above: The Hausdorff measure of the rational 
numbers £ [0, 1] is 0, whereas the Hausdorff measure of the real numbers e [0, 1] is 1. Therefore, 
the Hausdorff measure of a set is a proper measure. The Hausdorff measure of the fractal distri- 
bution in phase space M{t — > oo, to) is the same as that of M(t n ), W(E, N,Vi). Measured by the 
Hausdorff measure the phase space volume of the fractal distribution M(t — > oo,i ) is conserved 
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and Liouville's theorem applies. This would demand that thermodynamics could distinguish be- 
tween any point inside the fractal from any point outside of it independently how close it is. This, 
however, is impossible for any macroscopic theory which has only macroscopic information where 
all unobserved degrees of freedom are averaged over. That is the deep reason why the box-counting 
"measure" must be taken and is a further origin for irreversibility. 
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